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PREFACE 


The  study  reported  in  this  Memorandum  is  an  outgrowth  of  RAND's 
continuing  interest  in  the  electromagnetic  effects  of  nuclear  explosions . 
It  should  be  useful  to  students  of  the  theoretical  aspects  of  systems 
designed  to  function  in  the  environment  of  such  explosions. 
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SUMMARY 

The  close-in  electromagnetic  fields  produced  by  deflection  in 
the  earth's  field  of  Compton  electrons  from  a  nuclear  explosion  are 
analyzed.  Maxwell's  equations  in  spherical  coordinates  are  solved 
by  an  expansion  in  the  perturbation  fields,  taking  into  account  the 
space  and  time  dependence  of  the  conductivity  and  Compton  current. 

The  field  structure  is  determined,  and  it  is  shown  that  the  peak  change 
in  field  is  only  10  per  cent. 
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1 

I.  INTOOgJCTION 

It  is  well  known  that  nuclear  explosions  may  produce  appreciable 

(1  2) 

electromagnetic  signals.  '  Hie  theory  has  been  developed  for  the 
electromagnetic  radiation  from  nuclear  explosions  at  high  altitudes  or 
in  space ,  and  also  for  Compton  electron  interactions  with  the 
earth's  field. ^  However,  most  of  the  theory  deals  with  electric 
fields  or  dipole  radiation  fields.  This  Memorandum  considers  the  near 
magnetic  fields.  The  principal  mechanism  treated  is  Compton  electron 
interaction. 

The  approximations  vised  in  the  analysis  restrict  the  solution  to 
early  times.  However,  the  fields  may  be  expected  to  be  strongest 
shortly  after  the  arrival  of  the  gamma  rays,  so  these  results  represent 
the  most  significant  portions  of  the  field.  The  nature  of  the  analysis 
is  such  that  the  results  should  be  upper  bounds  on  the  actual  fields. 
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II.  CLOSE-IN  ELECTROMAGNETIC  FIELDS  PRODUCED  BY  COMPTON  ELECTRONS 

The  analysis  will  toe  devoted  primarily  to  the  possible  exclusion 
of  the  earth's  magnetic  field  toy  the  electrons  produced  toy  the  nuclear 
blast.  To  avoid  great  complications,  the  effect  of  the  proximity  of 
the  conducting  earth  will  toe  neglected.  If  the  earth  is  highly  con¬ 
ducting,  the  image  charges  induced  in  it  will  tend  to  cancel  the 
horizontal  electric  and  vertical  magnetic  fields  near  the  surface, 
and  to  augment  the  vertical  electric  and  horizontal  magnetic  fields 
toy  a  factor  near  2.  These  effects  may  toe  considered  in  a  future 
Memorandum. 

If  the  proximity  of  the  earth  is  neglected,  the  gamma  rays 
produced  by  the  blast  will  toe  emitted  primarily  in  the  radial  direction. 
These  gamma  rays  will  toe  scattered  and  absorbed  in  the  atmosphere, 
producing  Compton  electrons.  We  shall  only  consider  low-altitude 
blasts,  for  which  the  variation  of  atmospheric  density  with  altitude 
will  toe  neglected.  Since  the  scale  height  is  about  8.4  km,  and  we 
are  interested  in  heights  below  3  km,  this  is  a  reasonable  assumption. 
Under  these  circumstances,  the  atmospheric  conductivity  will  only  toe 
a  function  of  time  and  distance  from  the  blast.  In  MKS  units,  the 
conductivity  is  given  toy: 

o(t,r)  =  e  u)e  n(t,r)  (l) 

where  e  is  the  electron  charge  in  coulombs,  uue  the  electron  mobility 
in  (meters/sec) /(volts/meter),  and  n(t,r)  is  the  number  of  electrons 
per  cubic  meter  at  the  distance  r  from  the  burst  point  at  the  time  t. 
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The  proper  value  o f  mobility  to  employ  is  somewhat  uncertain. 

The  electrons  are  rapidly  slowed  down  by  collisions,  and  then  they 

become  attached  to  oxygen  via  a  three-body  reaction.  The  duration 

-8 

of  this  process  is  about  1  shake  (10*  sec).  Measurements  of  the 

electron  mobility  in  air'^”^"^  indicate  that  the  mobility  and  also 

the  electron  attachment  rate  are  energy- dependent. 

The  data  are  given  as  functions  of  the  ratio  of  electric  field 

strength,  E,  to  pressure,  p.  At  atmospheric  pressure,  and  the 

field  strengths  estimated  to  be  produced  by  a  1-MT  blast,  E/p  is  so 

low  that  there  simply  are  no  published  measurements  in  this  region, 

and  it  is  necessary  to  extrapolate  the  curve.  OTiere  is  appreciable 

curvature  in  the  mobility  curve  at  the  lowest  values  of  E/p  for  which 

there  are  data,  so  the  extrapolation  may  be  subject  to  error.  There 

is  also  some  evidence  of  a  cutoff  in  the  collision  process,  which 

causes  the  mean  energy  of  the  electrons  to  be  higher  and  changes  the 

effective  mobility.  In  all,  estimates  of  the  mobility  may  be  in  error 

by  a  factor  of  2  or  3.  After  discussions  with  R.  BJorklund  and 

W.  Karzas  on  the  data  extrapolation  process  an  average  value  of  .20 
2 

meters  /volt-sec  has  been  selected  for  the  calculations.  The  electron 
density  produced  by  a  1-MT  burst  has  been  calculated  from  Eq.  (4)  of 
Ref.  6j  and  the  conductivity  at  selected  values  of  distance  and  time 
is  listed  in  Table  1.  The  time  is  in  shakes,  measured  from  the 
initiation  time  r/c,  and  the  distance  is  in  kilometers. 

These  data  have  been  deduced  from  the  expression 

G(r)  g(t  -  £) 


a  = 


(2) 


k 


Tibia  1 

ATH5SPEERIC  CONDUCTIVITY  (MBO/METTO)  PRODUCE)  BY  A  1-MT  EXPLOSION 


Tine, 

Shakes 


t. 


Distance, 


r, 


fron  burst  point,  km 


0.5 

0.75 

1.25 

m 

1.75 

2.5 

3.0 

1 

2.9(-2) 

6. 3(-3) 

1.6(-3) 

4.0(-4) 

l.K-4) 

3>7(-5) 

6.o(-6) 

9.3(-7) 

2 

1.7 

3-3 

7-0 

1.9 

5.0 

1.4 

4.0 

6.7 

1.2(-6) 

5 

1.6 

2.9 

6.3 

1.7 

5.0 

1.5 

4.3 

6.0 

1.3 

10 

1.3 

2.4 

5-3 

1.4 

4.3 

1.4 

4.3 

5.0 

1.1 

20 

9.0(-2) 

1.7 

3.8 

1.1 

3.3 

1.1 

3.8 

4.5 

8.3(-7) 

50 

3-8 

7.7(-3) 

2.2 

6.7(-V) 

2.2 

8.0(-5) 

3.0 

4.0 

6.3 

100(1  u 

sec) 

2.3 

5.0 

1.4 

4.8 

1.7 

6.3 

2.5 

3.8 

5.7 

200 

1.6 

3-7 

1.1 

3.7 

1.3 

5.0 

2.1 

3.6 

5.0 

500 

6.7(-3) 

2.1 

1.6 

5.8(-5) 

2.4 

1.1 

2.2 

3.7 

1.0(3) 

1.6 

4.0(-4) 

fm 

4.0(-5) 

1.7 

7.8(-6) 

4.1(-6) 

1.3 

2.8 

2.0(3) 

9.7(-5) 

2.7(-5) 

9.3  (-6) 

M(-6) 

2.7(-6) 

1.9 

1.4 

7.3(-7) 

2.3 

5-0(3) 

5.2(-6) 

1. 9( -6) 

1.0(-6) 

7-3(-7) 

5.7(-7) 

5.0(-7) 

3.7 

1.7 

1.0(4) 

5.2 

1.9 

. 

6.7(-7) 

4.3 

5.1 

2.5 

1.8 

1.2 

Numbers  in  parentheses  denote  powers  of  10 
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G(r)  -  e-(r/X)/(r/x)2  (3) 

s(t  '  f>  "  «o  f  °  d  t#  f(t  -  |  -  t')  e*P  t#  (4) 

where  \  is  the  mean  path  for  removal  of  the  gamma  rays  (300  meters), 
p  is  the  attachment  rate  of  the  electrons  (10°  sec*  ),  gQ  is  a 
constant  which  is  related  to  the  peak  value  of  the  conductivity,  and 
f (t)  is  the  normalized  production  rate  of  the  gamma  rays.  For  the 
listed  data,  a  good  value  for  gQ  is  130,  while  the  time  dependence  of 
g  is  rather  complicated.  A  characteristic  time  to  represent  the  rise 
of  g  is  2  shakes,  while  the  decay  time  is  on  the  order  of  20  shakes. 

In  these  units,  the  Maxwell  equations  are: 


Vx  B  =»  H  +  h(cte‘  +  ‘j) 

(5) 

0  « 

(6) 

where  ii  is  the  permeability  (4n  x  10*^),  c  the  velocity  of  light,  J 
is  the  Compton  current  density  ( amps /meter  ),  E  the  electric  field 
(volts/meter)  and  B  the  magnetic  induction  (webers/meter  ). 

The  Compton  electrons  will  be  deflected  by  the  local  magnetic 
field.  Practically  all  of  the  effect  is  due  to  the  primaries,  since 
the  secondaries  undergo  very  little  deflection,  tod  the  deflection 
contributions  in  the  two  directions  perpendicular  to  the  magnetic 
field  roughly  balance.  Following  Karzas  and  Latter ' ^ ^  the  Compton 
current  is  given  by: 


J 


max  (t  -  £,  £) 

"j(t,r)  ~  JQ  G(r)  f  (t  -  |)  f  dt'v(t')C  (7) 

o 

where  the  velocity  of  the  Comptons  has  been  set  equal  to  c  in  the 
production  function.  JQ  is  a  constant,  G  and  f  are  as  in  Eq.  (2). 

The  upper  limit  of  the  integral  is  equal  to  the  larger  of  t  -  ^  and 
—  ,  where  R  is  the  range  of  the  Comptons.  Since  R  ~  1  meter  at  sea 
level,  the  approximation  R/\  «  1  has  been  made  (actually,  R/>  ~  .01). 

For  the  listed  data,  the  constant  JQ  is  6  x  10~\ 

The  Comptons  are  scattered  radially  from  their  point  of  origin. 

Take  a  system  of  spherical  coordinates,  with  the  origin  at  the  burst 
point,  and  the  polar  axis  along  the  earth's  magnetic  field;  then  the 
velocity  of  the  deflected  electrons  is  approximately: 

e  B  7  v  t' 

v  ~  yo  ir  ♦ - ~  “  sin  9  *0  (Q) 

Thus  the  deflection  is  in  the  azimuthal  direction,  has  a  latitude  dependence 
sin  0,  and  the  entire  system  is  azirauthally  symmetric. 

Under  these  conditions,  the  Maxwell  equations  split  into  two  groups, 
one  involving  the  radial  electric,  latitudinal  electric,  and  azimuthal 
magnetic  fields;  the  second  involving  the  azimuthal  electric,  latitudinal 
magnetic  and  radial  magnetic  fields.  The  first  is  driven  by  the  radial 
current,  the  second  by  the  azimuthal  current.  If  the  radial  current  is 
only  a  function  of  distance  and  time,  as  it  is  to  this  approximation, 
the  first  group  of  equations  reduces  to  Poisson’s  equation,  involves 
no  magnetic  fields,  and  will  not  be  considered  further. 

In  the  spherical  coordinates,  the  second  group  of  equations  may 


be  written  as: 
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,  >  3B 

- ,±_  4-  sin  e  Ew  =  -  rr-r- 

r  sin  %  ae  <p  at 


-  It  r  E„ 


-  ih 

r  a?  *'  y  ~  at 
5  b  ,  5l 


(9) 

(10) 


r  |rBe'  r  ST  =  1  rif  +  11  G(r)  [g(t  ‘  c)  Ej 

•1  < 


+  h(t  — )  sin  9 
c 


e  B  V  v 

h(t  -  f)  =  J0f(t  -  f) 


max 


(t 


(12) 


Our  task  is  now  to  solve  these  equations  with  the  appropriate  initial 

x* 

conditions.  For  t  <  -,  the  electric  field,  conductivity,  and  current 
vanish,  while  the  magnetic  field  is  uniform.  The  angular  dependence 
is  established  by  the  factor  sin  0  in  the  current.  A  representation 
which  integrates  the  first  two  equations  is: 


II 

sin  6  5  \ 

r  5t  F(t’r) 

(13) 

Be  = 

-B  sin  6  +  — ^  F(t,r) 

(14) 

B  = 
r 

B  cos  0  -  -  C°S-~  F(t,r) 
r 

(15) 

where  F  satisfies  the  wave  equation 

■US’ 

1 

II 

h  J  +  ^  G<r>  [s(t  - 1)  I; +  rh(t  -  f)  ] 

(16) 

with  the  initial  conditions 


0 


8 


(17) 


'  -I? 


r 

"c 


P  Is  related  to  the  vector  potential  function. 

Introduce  as  new  variables  the  distance  r, 
labeled  p,  and  the  delayed  time  t  =  t  -  The 


pressed  in  these  variables  is 


which  will  be  re¬ 
wave  equation  ex- 


b2?  2F  2  b2  F 

i?  ’  p2  ’  * 


U  Gr(p )  [g(T)  ^  +  P  h(t)J 


(18) 


with  the  initial  conditions: 


F 


c3F 

8t 


0 


T  =  0 


(19) 


The  solution  obtained  by  Karzas  and  Latter' ^  is  equivalent  to 
neglecting  the  left  side  of  this  equation.  We  wish  to  determine 
the  effect  of  this  neglect. 

For  this  purpose,  consider  characteristic  lengths  and  times. 

The  p  dependence  is  characterized  by  the  mean  free  path  \(300  meters). 
Use  current  and  conductivity  are  characterized  by  times  between  2  and 
20  shakes,  or  equivalent  lengths  6  and  60  meters.  Therefore,  the 
t  variation  is  fast  compared  to  the  p  variation.  Let  us  introduce 
dimensionless  variables  by: 


p  =  X  x 

X  =  300  meters, 

1  s  X  s  10 

(20) 

t  =  T  y 

T  =  10  sec 

0  <  y  s  2000 

(21) 

where  the  inequalities  indicate  the  expected  range  of  values  of 
significance  of  the  variables.  The  wave  equation  becomes: 


r  $ '  7)  ‘ 2  =  M  c  x  0(x)  [g(y)  % 

+  X  T  x  h(y)J  (22) 

The  essential  point  of  the  transformation  is  that  c  T/X  =  .01.  It 

may  be  expected  that  derivatives  with  respect  to  x  and  y  are  of  the 

same  order  of  magnitude.  If  ve  are  not  too  close  to  the  origin,  say 

-2 

x  >  1,  the  term  involving  x  is  on  the  same  order  of  magnitude  as 
the  x-derivative  term.  It  is  therefore  plausible  to  expand  the 
solution  in  powers  of  c  T/X,  and  keep  only  the  first  two  terms.  The 
equation  splits  as  follows: 

F  *  F0(x,y)  +  Hi  Fx(x,y)  (23) 


sra?  + 1  * cX  G<x)  «(y>  sr  a  - 1 


These  equations  may  be  solved  explicitly,  since  they  are  first- 
order  linear  equations  in  3F/dy.  Introduce  a  function  a(y)  by  the 


relation 
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a(y)  =  |  M  C  X  g(y)  =  188  X  g  (y)  (27) 

The  function  a(y)  represents  half  the  ratio  between  the  characteristic 
length  a  and  the  characteristic  length  l/p,  cg(y)  associated  with  the 
magnitude  of  the  conductivity.  For  large  conductivities,  a(y)  is  a 
large  number.  For  the  data  represented  in  Table  1,  a(y)  has  a 

c  z 

maximum  value  2  x  Hr  at  y  5  (5  shakes),  and  exceeds  Hr  out 
to  y  =  1000  (10  microseconds). 

The  solution  of  Eq.  (24)  subject  to  the  initial  condition  Eq. 

(26)  is:  x 

-a(y')  f  dx"  G(x" ) 

y  x  J  x  / 

FQ(x,y)  =  -  |  u  c  X2T  j  <3y/  h(y')  J  dx7  x'  G(x')  e 

(28) 

The  lower  limit  of  the  x1  integration  has  been  set  equal  to  zero  to 
prevent  the  solution  from  becoming  exponentially  large  at  x  =  0. 

The  lowest-order  field  components  are  given  by: 


V(x,y)  =  "  c  Beo(x’y)  = 


p  c  X 
2  x 


fc(y) 


x 

j  dx'  x'  ij(x')  e 


X 

f  dx" 
•C  / 


G(x"  ) 


(29) 


dx'  x'  G(x') 


dx"  G(x") 
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The  factor  h(y)/g(y)  is  the  Karzas-Latter  solution.  It  is  therefore 
necessary  to  study  the  properties  of  the  remaining  factor  of  Ec.  (30), 
which  will  he  called  W(x,y).  First,  since  a(y)  and  G(x)  are  positive, 
W(x,y)  is  positive.  Therefore,  the  Karzas-Latter  result  that  the 
magnetic  field  is  augmented  rather  than  excluded  remains  valid  to  this 
order  of  approximation.  Second,  an  integration  by  parts  brings  W 
into  the  form: 

x 

x  -a(y)  J  dx"  O(x') 

W(x,y)  -  1  -  j  J  dx'  e  *'  (30) 

o 

The  integral  appearing  here  is  positive,  so  W(x,y)  is  less  than 
unity.  Therefore,  the  Karzas-Latter  solution  is  always  greater  than 
the  solution  obtained  here. 

A  detailed  study  of  the  function  w(x,y)  is  presented  in  the 
Appendix.  A  simple  approximate  form,  correct  to  about  1  per  cent, 
has  been  obtained.  Define  a  characteristic  length  x^y)  as: 


a(y) 


*l(!  +  xx)  e 


*1 


(3D 


For  the  data  of  Table  1,  x^  is  between  4  and  8  over  the  range  1  shake 
to  10  microseconds.  The  length  x^  is  the  radius  of  the  ionization 
sphere,  which  is  thus  between  1200  and  2400  meters.  The  function  W 
is  closely  represented  by  _ _ 


W(x,y)  =  1 

^(y) 


X 


X  <  Xx(y) 
x  a  x1(y) 


(32) 
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so  the  Karzas-Iatter  result  is  quite  accurate. 

The  functions  x^(y)  and  h(y)/g(y)  are  plotted  in  Figs.  1  and  2 
for  the  data  of  Table  1,  and  the  time  range  y  <  7*  The  current- 
conductivity  ratio  has  been  normalized  by  dividing  by  cB,  producing 
a  dimensionless  number  which  can  be  regarded  as  a  field  susceptibility 
(perturbation  magnetic  field  produced  by  electrons  deflected  by  a  unit 

magnetic  field).  For  y  >  7,  the  ratio  Jq/c  B  aQ  has  the  constant  value 

* 

0.10.  This  constancy  is  established  by  the  exponential  character  of 
the  decay,  which  may  be  expected  for  most  nuclear  explosions.  This 
decrement  constant  is  a  very  gradual  function  of  the  yield.  The 
effective  radius  x^(y)  decreases  approximately  linearly  as  y  exceeds 
7.  The  slope  of  the  line  is  approximately  .05  until  values  of  x^ 
less  than  unity  are  reached,  which  for  the  listed  data,  takes  place  at 
y  »>  200. 

For  the  1-MT  explosion,  the  maximum  value  of  x^(y)  is  8.  Therefore, 
for  distances  greater  than  x  =  8  (2.4  km),  the  azimuthal  electric  and 
latitudinal  magnetic  fields  rise  to  the  maximum  value  8  Jq/oq  x,  and 
then  decrease  linearly  with  time  with  slope  .05.  For  x  less  than  8, 
the  fields  display  a  plateau,  which  begins  when  x^(y)  first  reaches  x, 
and  lasts  until  the  second  crossing. 

The  numerical  constant  which  determines  the  asymptotic  field 
strengths  is  given  by: 

* 

Karzas -Latter  use  0.05  for  this  number.  Discussions  with  R. 
Bjorklund  and  W.  Karzas  concerning  the  interpretation  of  the  experimental 
data  have  established  0.10  as  more  accurate.  Hence,  the  fact  that  our 
fields  are  3*3  times  Karzas-Latter 's  upper  bound  is  a  consequence  of  the 
choice  of  constants,  and  the  analyses  are  consistent. 
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J 

c  c  B 


(35) 


For  Compton  electrons,  y  =  2,  vq  ~  c.  Here  the  gamma  ray  spectrum  of 
the  device  has  been  employed.  The  mean  range  of  the  Compton  is 
taken  as  1  meter,  corresponding  to  .44  MEV.  This  number  is  obtained 
from  the  mean  energy  resulting  by  integrating  the  Klein-Nishina 
scattering  formula  over  the  primary  spectrum.  The  mobility  is  .20, 
and  q  is  the  ratio  of  the  number  of  secondary  to  primary  electrons 

I4. 

(3  x  10  ).  The  insensitivity  of  the  number  to  yield  is  apparent. 

(Dr.  R.  Bjorklund  provided  the  required  data  for  these  numbers . )  From 
the  numbers  presented  here,  it  follows  that  the  peak  change  in  the 
vertical  field  is  on  the  order  of  10  per  cent  of  the  earth's  magnetic 
field. 

The  field  line  structure  will  be  perturbed  by  the  current  flow. 
The  equation  for  the  lines  of  the  perturbation  field  is 


dr  _  _  2  cos  9  F 

rd©  "  ABg  ~  "  sin  0  r  dF/dr 


which  may  be  integrated  to  yield: 


(34) 


F(r,t)  sin2©  =  F(rg,t)  (35) 

where  rg  denotes  the  value  of  r  at  which  a  particular  field  line  crosses 
the  equatorial  plane  (9  =  £). 

The  function  F  must  be  obtained  by  numerical  integration.  Equation 
(3t>)  represents  the  field  lines  of  a  particular  instant  of  time,  where 
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F  has  been  specified  as  a  function  of  r  and  of  t  -  £  by  the  wave 
equation,  Eq.  f 18) ,  and  its  approximate  solution  Eq.  (28).  Accordingly, 
the  integration  is  quite  complicated.  The  results  are  presented  for 
a  time  t  *  10  microseconds  after  the  initiation  of  the  blast,  at  which 
time  the  leading  edge  is  at  3  km  (10  x  units).  The  function  F  should 
not  be  taken  too  seriously  for  small  values  of  x,  although  the  general 
behavior  is  correct.  Figure  3  is  a  plot  of  F(x)  at  ID  microseconds. 

Since  F  has  a  maximum,  a  field  line  must  cut  the  equatorial 
plane  twice.  Therefore,  the  lines  close  on  themselves.  The  magni¬ 
tude  of  the  field  strength  varies  along  a  given  line,  being  largest 
at  the  outer  crossing.  The  perturbation  field  lines  are  plotted  in 
Fig.  4.  These  lines  are  arranged  to  have  equal  increments  of  field 
strength  at  the  outer  crossing. 

The  lines  are  strongly  crowded  near  the  front,  and  spread  as 
we  move  inward.  The  line  drawn  closest  to  the  front  is  actually  that 
on  which  the  perturbation  field  strength  at  the  outer  crossing  has 
its  largest  value,  10  per  cent  of  the  earth's  field.  There  are  as 
many  lines  between  that  line  and  the  front  as  there  are  within  that 
line.  Near  the  front,  the  field  is  predominantly  in  the  latitudinal 
direction,  while  at  long  distances  within  the  front  the  field  becomes 
nearly  radial,  except  at  the  equatorial  plane.  This  general  structure 
holds  for  all  times,  though  the  details  will  vary  as  the  front  expands. 

Next,  we  consider  the  higher  approximations  to  the  solution. 

Equation  (25)  may  be  solved  in  exactly  the  same  manner  as  Eq.  (24), 
and  the  same  type  of  approximations  made.  The  result  is  that  the 
field  near  the  front  is  affected  only  very  slightly  by  the  second  approxi¬ 
mation.  Farther  back  the  effect  is  stronger,  but  the  field  is  considerably 
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Flg'4  The  Per,urbotion  field 


lines  ot  t=io  jisec 
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weaker.  The  higher  approximations  will  contribute  only  about  1  per 
cent  of  the  total  field,  which  ie  much  less  than  the  uncertainties  in 
the  constants  and  in  the  validity  of  the  general  model.  This  result 
means  that  the  perturbation  analysis  is  self-consistent.  The  peak 
field  is  effectively  independent  of  yield,  hut  the  "equivalent  fire¬ 
ball  radius"  depends  logarithmically  on  yield. 

In  summary,  the  close-in  electromagnetic  fields  produced  by 
deflection  of  Compton  electrons  in  the  earth's  field  have  been  analyzed, 
and  it  has  been  shown  that  there  is  no  significant  exclusion  or  augmenta¬ 
tion  of  the  earth's  field. 
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APPENDIX 

The  function  W(x,y)  has  been  defined  by  Eqs.  (3)  and  (29)  as: 


W(x,y) 


x 

J  dx'  x#  G(x') 
0 


dx"  G(x* ) 


(36) 


G(x)  =  e"x/x2 


(37) 


The  argument  of  a  may  be  dropped,  since  W  is  a  function  of  x  and  of  a 
only.  As  x'  approaches  zero,  the  integral  in  the  exponent  tends  to 
positive  infinity,  and  the  integrand  tends  to  zero.  Bie  integral  may 
be  rewritten  as: 


W(x,y) 


X 

r  x  1 

I  dx'  exp  - 
J 

x'  +  log  x'  +  a  /  G(x')  dx' 

0 

V 

(38) 


The  expression  in  brackets  is  large  for  x'  very  large  or  very 
small.  It  therefore  possesses  a  minimum,  which  may  be  called  x^. 
Call  the  bracketed  expression  t(x' ).  Differentiating  and  setting 
the  result  equal  to  zero  yields: 


t'txj)  -  1  +  “■  -  a  (Kx^  -  0 

X1 

a  -  x1(l  +  x^  e 


(39) 


(40) 


SI 


Equation  (40)  gives  a  as  an  explicit  function  of  3^,  or  as  an 
implicit  function  of  a.  The  solution  of  this  implicit  relationship 
has  been  plotted  in  Fig.  1  against  y,  using  the  relation  between  a 
and  y  defined  in  the  text.  The  integral  in  the  exponent  may  be  written 
as: 


dx' 


J  G(x* )  dx'  -J  G(x')  dx' 

X '  X 


=  I(x')  -  I(x) 


(41) 


Figure  5  gives  I  versus  x  and  a  versus  x^ .  For  a  positive, 
Eq.  (40)  has  only  one  solution  with  x^  positive. 

The  value  of  W  will  depend  strongly  on  whether  or  not  the 
minimum  point  x^^  is  within  the  range  of  integration.  If  it  is 
within,  the  Integrand  takes  its  largest  value  at  x  =  x^,  and  decreases 
rapidly  as  x  moves  away  from  x^  in  either  direction.  In  this  case, 
the  integrand  may  be  approximated  by  the  well  known  method  of  steepest 
descent.  3he  exponent  is  replaced  by  the  first  three  terms  of  its 
Taylor  series  expansion  about  x  =  x^,  and  the  resulting  expression 
yields  a  close  approximation  to  the  integral.  There  results: 


W(x,y) 


1  ♦'  (Xj)  (*'  -  x^ 
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where  Jr*  denotes  the  second  derivative  of  evaluated  at  x. ,  and 
the  error  function  erf  is  defined  hy: 


erf  z  =  — —  f  du  e“u  (43) 

Equation  (39)  yields  for  the  second  derivative: 

♦  '  (xx)  «  1  +  J-  +  Jg  (44) 

1  *1 

For  greater  than  2,  the  argument  of  the  second  error  function  ex¬ 
ceeds  2.34,  and  the  function  is  between  .999  and  1.  The  departure 
from  unity  may  be  neglected. 

The  first  error  function  is  zero  at  x  ■  x^,  and  increases  rapidly 
to  unity  as  x  increases.  Except  for  a  small  shoulder  immediately 
above  x^,  the  sum  of  error  functions  may  be  set  equal  to  2.  The 
remaining  factors  have  been  studied  numerically.  The  variation  of 
\|r  with  x^^  arising  from  th.  term  l(x)  tends  to  cancel  the  variation 
of  the  error  function  terms  near  x  ■  x^.  For  x  significantly  larger 
than  x^,  l(x)  becomes  negligible.  The  other  parts  of  W  depend  on 
x^  only.  The  details  yield: 


2k 


x  >  x1(y) 


(45) 


For  x  less  than  x^,  the  minimum  is  not  vithin  the  range  of 
integration.  The  largest  values  of  the  integrand  now  occur  at  the 
end  point  x'  =  x.  Again  expanding  the  exponent,  but  now  around  the 
end  point,  there  results: 


W(x,y) 


t'(x)(x'-x)  +|  *'(x)(x'-x)2 


(46) 


1 

5 


J± 

♦'(x) 


>/2p00 


imiU 

♦'(*)  /  / 


Jim ' 

•f  2)*  (x) 


(47) 


Here,  the  arguments  of  both  error  functions  are  large.  The  difference, 
multiplied  by  the  exponential,  yields  an  algebraic  function.  A 
numerical  study  shows  that  the  resulting  expression  varies  by  only 
a  few  per  cent  for  x  ranging  from  0  to  7,  which  is  the  largest 
value  of  x^  in  the  detailed  problem.  Thus,  the  expression  may  be 
set  equal  to  unity,  since  it  must  be  continuous.  Therefore,  the 
function  W  is  given  by: 
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W(x,y)  -  1 

^(y) 

“  x 


x  <  x^y) 
x  >  x1(y) 


m 


This  expression  is  in  error  by  only  a  few  per  cent  over  the  complete 


range 
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